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SstlMitlng  the  Parameters  of  a Truncated  Gumma  Distribution  y 
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Unl-rersitjr  of  Washington 

1*  £aH£X»  A table  is  giren  to  slniplify  the  estimation  of  the  parametexti 
of  an  incomplete  gsoma  or  Type  III  distributioiu  A now  procedure  is 
also  suggested  to  estimate  the  parameters  of  a truncated  guna  distri* 
bution*  This  method  Is  also  applioable  to  a nuaber  of  other  truncated 
distributions^  whether  the  truncation  is  in  the  tails  or  the  center  of 
the  distributico. 

2.  Introduction.  Sereral  examples  hawe  been  giTon  recently,  tmploying  the 
incomplete  gmmui  or  Type  III  distribution  in  fitting  rainfall  data, 
e.g.  In  an  animal  population  study  It  was  found  that  the 

migration  pattern  could  be  fitted  by  this  type  of  distribution*  Fre> 
quently  in  such  aigration  studies  the  data  will  be  truncated,  i.e. 
obseanratians  will  begin  after  nigraticn  has  begun  or  oonolude  before  it 
has  stopped. 

The  parameters  of  the  gassm  distribution  are  often  estimated  by 
the  method  of  momants  in  such  cases  [see  for  example  [a]  pp.  121,  I25J 
despite  the  fact  that  Fisher  [5  3 showed  the  method  to  be  inefficient* 

To  fhcilitate  solution  of  the  maxlav  likelihood  equations  for  estima- 
tion of  the  parameters  in  the  untrunoated  case,  a simple  table  is  giwKU 
The  estimation  of  the  parameters  of  a truncated  gamma  distribution, 
hy  the  method  of  moments,  has  been  studied  by  C<^isn  [hj.  Since  the 
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lnt«gr«l  of  the  probability  doneity  cannot  bo  expreosod  in  closed  fozia^ 
OYon  the  mcnent  estimates  are  tedious  to  obtain;  no  attempt  has  been 
made  to  evaluate  their  Tarianeea  or  study  their  efficiencies.  After 
this  paper  vas  completed,  a nev  study  of  the  problem  mas  published  by 
Des  Ba4  [l3j«  He  glTSS  the  maxlaum  likelihood  equatione  for  a nuaber 
of  oases  of  truncated  and  censored  sam;d.es,  mainly  howavmr  under  the 
assui^tioo  that  the  third  standard  moment  is  known.  These  equations 
can  be  solved  only  by  iterative  methods.  In  this  paper  a new  method  of 
estimatian  of  these  parameters  is  Introduced  which  is  easier  to  apply. 
The  asymptotic  variance-covariance  matrix  of  the  estimates  is  deterainsd. 

3*  with  origin  knomn.  The  density  function  of  the  gMsma  dis- 

tribution may  be  written  in  the  form 

f ) 7^  * ^ * 

= 0 X < e 

The  parameters  are  frequently  transformed  so  that  the  distribution  is 
eiqpressed  as  a function  of  the  nwan,  variance  and  skownase.  Since  the 
corresponding  sample  quantities  do  not  efficiently  estimate  the  para- 
meters such  a transformation  appears  to  be  misleading. 

The  maTianni  likelihood  equations  have  been  given  b7  Fisher  [5]p  vis. 


W J - (i  - o)  = 0 

OS  a 


Sine©  the  parameter  c determinee  ths  region  of  poaitlve  demity, 

Uxat  equation  (3)  gives  the  maximum  llkollhood  estimate  of  c musi  be 

justified  in  a silently  different  manner  than  by  routine  csloulua.  If 

b > 1 this  is  easily  done;  if,  hovever,  b^  1 equation  (3)  does  net 

glre  the  raax1,imm  likelihood  of  c»  In  this  case  f(z)  is  monotone  dr« 

creasing  for  z ^ c and  s * nln  z.  is  the  maxi  min  likelihood  estimate 

1 ^ 

of  o« 


Ve  consider  first  the  case  vhere  the  origin  is  knownt  so  that  c 
nay  be  set  equal  to  serO|  without  loss  of  genereJllty  and  equation  (3) 
drops  out.  Letting 


In  z^ 


(1)  and  (2)  yield 

(4)  y-  (b)  » In  b - * In  z - 

81n«  . th. 

W,  It  is  easy  to  c<»struct  a table  of  y Cb)  and  solve  for  b b}'^ 
inverse  interpolation,  k enwll  tabulation  of  y (b)  is  given  in  Table  I; 
e more  oomplete  tabulation  is  available  In  mimeographed  form  from  ^he 
Laboratory  of  Statistical  Research,  Oniverslty  of  Vashlngton.  Thtu*e 
y(b)  and  its  first  and  second  differences  are  tabulated  for  b**I  («01)5f 
5(0.1)20}  20(1)100.  The  table  was  checked  by  sumnlng  colunna  in  thi« 
basic  tables  and  should  be  correct  to  one  figure  in  the  fifth  deeliaali. 


4*  letlaatlop  in  the  truncated  case  with  known  origin.  The  denv  ,ty  functlcm 
is  now  written 
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(4)  f (x)  = 0 r 

ss  0 9lawhtT0 
t 

vhsre  K(a,b)  » C •“**  x^^  dx 

i 

Th«  maxlaaim  Itkalihood  functions  nov  In^olTs  dsrliimtlTss  of  K vl-^ 
rsspset  to  s and  b rospsctlTsly;  a double  entry  table  would  be 
necessary  to  obtain  the  aaxlnua  likelihood  estlaates  of  a and  b and 
eTsn  this  would  InTolxe  double  inwerse  interpolation* 

In  lieu  of  this  another  netbod  of  estijMition  is  proposed.  Let 
the  n observations  be  grouped  by  desses  ^1^^)* 

(i=s  l,2,...r)  where  4^  “ ^r"^  denote  by 

the  nuaber  of  observations  falling  in  class  i»  i.e.  between 


^-1 


b-1  , ^i  / 

x*^  dx  i Ke  4^  (2b^) 


(6)  Inp^-lnp  =a(4^^^-4t)^(»»-l)  In-r^f-ln^ 

9i^i  1^1 

(is*  1^2,. ..r  - 1)  to  the  degree  of  approadjaation  indicated  by  (5). 

The  font  of  equation  (6)  suggests  estlBating  a and  b by  a 

least  squares  procedure,  with  q^  replacing  p^.  This  can  be  Justified 

as  an  approxiaate  prooedure  hy  the  following  results.  To  terns  of  order  ^ 

n 

(7)  £(ta  q^)  - In  - i 
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(S) 

(9) 


£ (In  - In  p^) 


2 ^ l-P, 


% 


(in  —I)  (in  — i)  sa  •»  “ 


ThM*  rMulte  can  ba  obtained  by  expanding  la(q^/p^)  «ln(if‘ — ^ — ) 

in  a Taylor  aarlee  (aaeualng  that  for  large  n»  0)  and 

Pr(q^  > 2p^)  nay  be  neglected).  It  la  alao  neceastry  to  uae  the  fict 


that  £ (q^  f a reault  eaally  obtainable  rrcai  the  veU-imotm 

reourrenoe  foraula  for  the  central  noaenta  of  the  Bultlnoadal  diatribe* 
tlon,  els. 


/^r-M 


Froa  thla«  In  fact.  It  aay  be  Inferred  that 

£(,^) 


Moreover,  the  Halting  dlatributlon  of  the  In  la  eaally  obtained 
frca  the  following  laaaa. 


lank  j ” l»2,...»r)  be  a aequanoe  of  racdoa  varlablea 

and  a^(l~l,2,...,r)  be  conatanta  auch  that  the  Joint  dlatrlbup* 


tlon  of 


^1 


(1*  l,2,...,r) 


tenda  to  the  lialtlng  dlatributlon  F(y^,y2>«**»yp  m n^oa. 


and  l»t  f(x)  be  of  class  in  the  neighborhood  of  /^t • » • 

vrt,thXi 


th«>  Zin)  . tJJ^  

cr-i  . Ti 


have  the  same  Joint  limiting  distribution. 

This  is  a consequence  of  the  general  theorems  on  stochastic  limit 
relationships  prored  by  Mann  and  Vald  Cs3  (cf.  their  Theorems  3 and  5-a 
trirlal  modification  is  howerer  required  since  our  f(x)  is  a function 
of  a single  real  rarlable  while  their  corresponding  g(x)  is  a function 
of  a yector-Talued  randoa  yarlable). 

Finally,  writing 

(10)  ■ In  - In  (l  »l,2,...,r  - 1) 

it  follows  that  the  yj^  are  asymptotically  multinormal  with  means 


(^i+1  • ^l) 


(b-1)  In  t In 


I 


i>l 


and  moment  matrix 


“Vii  Pj/  “VP2/ 

■) 


- 

"lf-2 


$ - 


“U/ 


0 


0 


Asymptotically  effloient  estiiaators  of  a and  b aro  foimd  by 


minimi  ling  the  quadratic  foxn 


> -1  ^ 

7 - f (y)*  7 - £(y) 

where  the  vector  ^=(y,y,...y  ) 

X 2 T^\. 


Since  the  true  valxiae  of  the  are  not  known*  it  le  neoeesary 
to  replace  the  p^  in  77^  by  their  eatimatee*  the  q^»  Introducing 
the  notation 
(11) 

(12) 


w a*  y — la  h ■/"  In  h 
1 ^1  1 . lt-1 


(13) 


^ IH  ■ ^1 
v^=  la  - 1*^ 


i-f-1 


the  eq;uatiati8  for  a and  b are 
(U) 


(«)  ^ {z  z^.i  v,}=|:  V4 

denoting  the  eXeaente  of  with  p*e  replaced  by  q*e) 


The  aoluticne  of  these  are 

(16)  ^ [(^  7?^w)(u*  w)  - (u*  ^ 

(17)  b =.  i [(u»  - (t  t)(u*  77j^^v) 

where  A = (u»  u)(^;^’^t)  - 

and  the  covariance  matrix  of  (a*  is 
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2 


^ (▼*  7?]  v)  - ^ (u*  7Tl  t) 
.-1 


.-1 


(u» 


Th«  estinatas  m and  b ara  found  by  dlract  alapla  routlna 
oalculations  axoapt  for  tha  datazslnation  of  7/^^  trm  7^.  This  nay 
ba  a tadious  prooass  imlaos  r la  lan  ■ Uowarar,  if  all  ara 
aqual  to  ^ » than 


2 

n 


This  is  aasily  Tsrifiad  by  diraet  nultiplioatian.  This  suggasts  that 
vhara  possibla  tha  ^ ^ should  ba  ohosan  so  that  tha  (axid  thus  ap- 
prcudnatalyy  tha  p^)  ara  acpial*  Tha  darioa  is  analogous  to  that  sug- 
gastad  by  Qunbal  £?J  and  by  Mann  and  ¥ald  in  applying  tha  X ^ 


"goodnass-of-fit*  tast. 

If  this  is  not  possibla,  lass  affloiant  astlaatas  oan  ba  obtainad 
by  utilising  only  tha  odd  (or  aran)  v^*s.  Tha  odd  v^'s  ara  sutually 
indapsndant  anong  thansalras  and  oonsaqusntly  7?^  and  raduea  to 
diagonal  natricasa 


$•  wWtwMH  ftyjgin-  If  tha  paranatar  o,  tha  origin,  is  un- 

known, than  tha  astination  problna  is  nora  difficult  vhathar  or  not  tha 

t 

A 

« 

* a»  f • 
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distribution  is  truncated.  Itfimtive  aothods  are  of  course  possiblo  in 
solving  (1)  (2)  and  (3)  vith  the  aid  of  Table  1,  i.Ov  for  the  untruncated 
case..  In  the  truncated  case  this  method  is  too  tedicis  to  have  much 
practical  value,. 

If,  in  the  tzunoated  case,  there  is  available  supp.\«nantaj;y  infoma- 
tion  so  that  the  restriction  0 < c < 4::  4^  utilised,  than  a 

procedure  similar  to  that  outlined  above  may  be  foUovec.  In  this  case 

(18)  In  - in  .( ^ -^(b-l)  In  — +ln  5^ 

(1  -l,2,...iwl) 

again  to  the  degree  of  approxlinatlon  indicated  by  (3).  With  the  restrlo> 
tion  noted  above,  it  is  adequate  to  vrlto 

(19)  In  p^  - In  •(  - 4^)  1®  f ^ ■ +(b-l)  o 

h 

1+1 

Defining  *0  above  least  squares  estimates  of  a,  b and  o 

may  be  found  in  an  exactly  analogous  procedure  to  that  of  Section  4. 

6.  Conclusion.  The  method  used  to  estimate  the  parameters  a and  b in 
Section  4 nay  also  be  applied  if  the  sample  is  drawn  from  a doubly 
truncated  gamma  distribution,  from  a singly  or  doubly  truncated  noxnal 
distribution  or  from  a beta  distribution  with  known  range,  either  trunoatsd 
or  not.  Methods  of  obtaining  the  maximum  likelihood  estimates  of  the 
parameters  of  a truncated  normal  distribution  are,  of  course,  well,  known. 
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and  axtansivo  tabniljitloiia  hare  bean  made  to  facilitate  ':he  determination 
of  auch  solutloDS  (e.g.,  compare  particularly  Bald  /llj)* 

The  method  outlined  abore  vould  also  be  useful  In  estimating  the 
parameters  of  the  normal  curre  where  there  are  systematic  gaps  in  the 
obserratloxis.  This  may  occur  particularly  In  time  distributions — an 
example  may  be  found  In  /l2ja  For  distributions  with  finite  but  unknown 
range,  howerer,  the  method  does  not  appear  to  be  satlsfactoryo 
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TABLE  I 


b 

1.0 

1.1 

1.2 

1.3 

1.1 

1.5 

1.6 

1.7 

1.5 

1.9 

2.0 

2.1 
2.2 

2.3 
2.A 

2.5 

2.6 

2.7 

2.8 

2.9 

3.0 

3.1 

3.2 

3.3 


y"(b) 

.57722 

.51906 

.47136 

.43156 

.39786 

.36898 

.34396 

.32208 

.30280 

.28567 

.27036 

.25660 

.24a6 

.23287 

.22257 

.21313 

.20446 

.19647 

.18907 

.18221 

.17583 

•16988 

.16431 

.15910 


,05816 

.04770 

.03980 

.03370 

.02888 

.02502 

.02188 

.01928 

.01713 

.01531 

.01376 

.012U 

.01129 

.01030 

.00944 

.00867 

.00799 

.00740 

.00686 

.00638 

.00595 

.00557 

•00521 

.00489 


L 


3o4 

.15421 

.00460 

3.5 

.14961 

.00434 

3.6 

.U527 

.00409 

3.7 

.14118 

.00387 

3.8 

.13731 

.00366 

3.9 

.13365 

.00347 

4.0 

.13018 

.00330 

4.1 

.12688 

.00314 

4.2 

.12374 

.00298 

4.3 

.12076 

.00284 

4.4 

.11792 

.00271 

4.5 

.11521 

.00259 

4.6 

.11262 

.00248 

4.7 

.uou 

.00237 

4.8 

.10777 

,00227 

4.9 

.10550 

.00218 

5.0 

.10332 

.009665 

5.5 

.093655 

.008013 

6.0 

.085642 

.006751 

6.5 

.078891 

.005765 

7.0 

.073126 

•004980 

7.5 

.068146 

•004346 

8.0 

.063800 

•003825 

8,5 

.059975 

.003392 

9.0 

.056583 

.003029 

9.5 

.053554 

,002722 


10.0 

.050832 

n.o 

.046143 

12.0 

.042245 

13.0 

.038954 

U.o 

.036139 

15.0 

.033704 

16.0 

.031575 

17.0 

.029700 

18.0 

.028035 

19.0 

.026547 

20.0 

.025203 

.C04689 

.003898 

.003291 

.002815 

.002435 

.002129 

.001875 

.001665 

.001488 

,001339 


